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where S q (u>) is the PDM strength function 



Report by N. Dinh Dang, V. Kim Au, T. Suzuki, and A. Arima 
on the Comment by V. Yu. Ponomarev* 

(October 24, 2001) 

In our recent work H , we calculated the El resonances in neutron-rich oxygen and calcium isotopes 
within a quasiparticle representation of the phonon damping model (PDM) in its PDM-1 version ^|J3) 
including the superfluid pairing interaction. It is claimed in the Comment p| that "the physical 
content of the PDM calculations is very doubtful" and the description of the pygmy dipole resonance 
(PDR) is "not justified" on "the quantitative level" . However, as seen below, the arguments, which 
the author of the Comment presents to prove his case, are cither wrong or irrelevant. 

The same Appendix 2D of Ref. [2] in the Comment ji[ states that the assumption of the 
equal particle-phonon coupling strength = f± "may be employed if the width of the strength 
function is small compared to the characteristic energies associated with systematic variation in 
the coupling matrix elements." This condition is satisfied in the region of El resonances con- 
■ sidered in [jl). However, this does not means that the matrix elements of the interaction part 

. V qisi = (qi \ J^ss'q Fss' a l a s'(Qq + Qq)\ s i) are the same, as claimed in the Comment 

Equation (1) of the Comment 0] is not the second moment for the phonon distribution within 
the PDM. The A:-th moment for the phonon distribution within the PDM is calculated (See, e.g., 
|^ ; Eq. (2.21) of §) as 

m ( q k) = S q (u)u> k ckj , fc = l,2,... (1) 

O 

9 ■ W = , ^ „ v (2) 

O 

The energy u> of the giant dipole resonance (GDR) is found as the solution of Eq. (2.39) in ||: 

Q-w q -P q (w) = 0, (3) 

where uj q is the unperturbed phonon energy (before the p/i-phonon coupling is switched on), and 
P q (u) is the polarization operator. The damping 7 9 (w) is calculated microscopically within PDM as 
the imaginary part of the analytic continuation of P q (E) into the complex energy plane E = u> ± ie. 
Its explicit expression within PDM-1 is given by Eq. (5) of (pairing not included) or Eq. (15) 

of |l| (pairing included). There is no way to equalize from (|l|) (k = 2) with Eq. (1) of the 
Comment Q]. Therefore, all discussions using Eq. (1) of the Comment Q] with the aim of attaching 
it to the PDM are irrelevant. 

The strength function (0) is not a Breit-Wigner (BW) distribution because the damping 7 9 (w) 
depends on the energy lo. Such form has been derived, for the first time, in Q using the analytic 
properties of the double-time Green function independently of any assumption on the coupling 
matrix elements. Consequently, the photoabsorption cross section of GDR within the PDM is not a 
Lorentzian either. The claim in ||] that a Breit-Wigner (Lorentzian) form is assumed or an ad hoc 
input for the strength function (photoabsorption cross section) of GDR within the PDM is simply 
wrong. 

One of the crucial features of the PDM is the use of realistic single-particle energies to construct 
the ph configurations (at zero temperature) together with the pp and hh configurations (at nonzero 
temperature) to which the GDR is coupled. Therefore a replacement of the realistic single-particle 
spectra with any other ones, such as the random values of E s used in the Comment no longer 
corresponds to the PDM. So the attempt in Q to imitate the results of the PDM using random 
values of E s , and an GDR energy Eq, which is not defined from (^|), is incorrect. 



*The Comment has been later rejected by PRC on January 31, 2002 
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It is by no mean obvious that the coupling constant should increase as the configuration space gets 
larger in the same nucleus. How the coupling changes is a problem to be discussed microscopically. 
Within the PDM, an increase of the space of ph pairs leads to a decrease of the parameter /i to 
preserve the same value for the GDR width. 

The aim of |jj is to use for the calculations of El resonances in neutron-rich isotopes the same 
set of two parameters (ui q , /i), whose values are chosen to reproduce the GDR in the corresponding 
double closed-shell nuclei. Therefore, the El resonances in the chains 16 ~ 24 0, 40_46 Ca, and 48_60 Ca 
were calculated using the parameters chosen for 16 0, 40 Ca, and 48 Ca, respectively. There is no 
reason why the values of /i for 40 Ca and 48 Ca should be the same. The results for GDR in 16 
have been obtained already within the enlarged space with fx = 0.6982 MeV. This value is kept 
unchanged throughout the chain of oxygen isotopes as has been mentioned above and in M. No 
change of the parameter occurs between 16 and 18 as incorrectly stated in the Comment 

The PDM-1 with its two phenomenologically selected parameters allows the comparison with the 
experimental data for only the average characteristics of the El resonances, such as the overall shape 
of the cross section, width, energy, and energy-weighted sum (EWS) of strength. It cannot describe 
such fine structure as the individual low- lying El states measured in 0. The EWS of El strength 
below 5 MeV for 40 Ca and 48 Ca are around 0.25 % and 0.52 % of the Thomas-Reich-Kuhn sum rule 
(TRK), respectively. They should be compared to the experimental values of (0.025±0.004) % of 
TRK for 40 Ca, and (0.29±0.04) % of TRK for 48 Ca reported in J|. The enhancement of strength 
at low energies in doubly closed-shell nuclei due to the spreading of GDR has been discussed in 
M. There is no such setting of B(E1) to zero for PDR within the PDM as claimed in [Q. The 
problem of double counting arises only when the structure of phonon is calculated microscopically 
within the random-phase approximation. The use of the structureless phonon and the parameters, 
which are selected so that the calculated GDR energy reproduces its experimental value, excludes 
any possibility for double counting within the PDM. 

Pairing is not included in the calculations of PDR in Ref. [5] of The aim of Rcf. [5] of Q 
is to see if the PDM is able to predict the existence of the PDR in neutron-rich nuclei, but not to 
reproduce the experimental data. Because of the absence of pairing, the parameter fx in Ref. [5] 
of Q was increased significantly in the region near the Fermi surface in neutron-rich isotopes. It is 
natural that such an increase overestimates the EWS of El strength in this region. This is not the 
case in the present work Q], where pairing is included, and the parameters of the model have been 
chosen to reproduce the GDR in double closed-shell nuclei. 

In conclusion, none of the statements of the Comment |Q is relevant. The discussions in the 
Comment Q are fruitless, and its conclusions are false. The Comment should not be published. 
If, nevertheless, the Editor will accept this Comment for publication, we would like to submit this 
report as our Reply to be published together with the Comment. 
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